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ABSTRACT 

O ' 

<n : 

P5 ' We study thermal-gravitational instability in simplified models for protogalac- 

iJ^ ■ tic halos using three-dimensional hydrodynamic simulations. The simulations 

00 ' started with isothermal density perturbations of various power spectra, and fol- 

lowed the evolution of gas with radiative cooling down to T = 10 4 K, back- 
ground heating, and self-gravity for up to ~ 20 cooling times. Then cooled 
and condensed clouds were identified and their physical properties were exam- 
ined in detail. In our models, the cooling time scale is several times shorter 
q ■ than the gravitational time scale. Hence, during early stage clouds start to form 

around initial density peaks by thermal instability. Small clouds appear first and 
they are pressure-bound. Subsequently, the clouds grow through compression by 
P H i the background pressure as well as gravitational infall. During late stage cloud- 

cloud collisions become important, and clouds grow mostly through gravitational 
merging. Gravitationally bound clouds with mass M c > 6 x 10 6 M are found 
in the late stage. They are approximately in virial equilibrium and have radius 
R c ~ 150 — 200 pc. Those clouds have gained angular momentum through tidal 
^ ■ torque as well as merging, so they have large angular momentum with the spin 

parameter (A s ) ~ 0.3. The clouds formed in a denser background tend to have 
smaller spin parameters, since the self-gravity, compared to the radiative cooling, 
is relatively less important at higher density. The H2 cooling below T = 10 4 K 
does not drastically change the evolution and properties of clouds, since it is much 
less efficient than the H Lya cooling. The slope of initial density power spectrum 
affects the morphology of cloud distribution, but the properties of individual 
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clouds do not sensitively depend on it. We point limitations of our study and 
mention briefly the implications of our results on the formation of protoglobular 
cluster clouds in protogalactic halos. 

Subject headings: gravitation — hydrodynamics — instabilities 

1. Introduction 

Thermal instability (TI) is one of key physical processes in astrophysical environments 
where optically thin gas cools radiatively and condenses (Field 1965). It has been applied to 
explain, for instance, the multiple phases of interstellar gas (e.g., Field, Goldsmith & Habing 
1969; McKee & Ostriker 1977), the formation of globular clusters (e.g., Fall & Rees 1985), 
cooling flows in clusters of galaxies (e.g., Nulsen 1986), and the generation of turbulent flows 
in the interstellar medium (ISM) (e.g., Koyama & Inutsuka 2002; Kritsuk & Norman 2002). 
Cold dense clumps or clouds that are confined by the background pressure can form in a hot, 
radiatively cooling medium via TI (e.g., see Burkert & Lin 2000). In the simplistic picture of 
TI, those overdense regions undergo a quasi-static compression in near pressure equilibrium 
(Field 1965). However, this isobaric condensation occurs only when the clouds are small 
enough to adjust to pressure changes faster than the gas cools. According to numerical 
simulations of the collapse of thermally unstable clouds (e.g., David, Bregman & Seab 1988; 
Brinkman, Massaglia & Miiller 1990; Malagoli, Rosner & Fryxell 1990; Kang, Lake & Ryu 
2000; Baek, Kang & Ryu 2003), the clouds may undergo a supersonic compression when the 
cloud size is comparable to the cooling scale, Z coo i (the distance over which a sound wave 
travels in a cooling time), while the clouds much larger than this scale cool isochorically. 

It has been suggested that the TI could be responsible for the formation of protoglobular 
cluster clouds (PGCCs) in protogalactic halos, which can explain the origin of old halo 
globular clusters (e.g., Fall & Rees 1985; Kang, Lake & Ryu 2000). Among many models 
of globular cluster formation (e.g., see Parmentier et al. 1999; Kravtsov & Gnedin 2005), 
this model based on TI is classified as a secondary model in which the two-phase of cold 
dense clouds in hot background gas was developed through TI in protogalactic halos and 
the condensed clouds further collapsed to become globular clusters. We previously studied 
the development of TI in detail using one and two-dimensional numerical simulations with 
spherically symmetric and axisymmetric isolated gas clouds in static environments of uniform 
density (Kang, Lake & Ryu 2000; Baek, Kang & Ryu 2003). However, according to the 
current paradigm of cold dark matter models of structure formation, large protogalaxies 
comparable to the Milky Way formed via hierarchical clustering of smaller systems. Hence, 
inevitably density perturbations should exist on a wide range of length scales inside the 
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protogalactic halos, and an ensemble of clumps emerge. In addition, although the thermal 
process initiates the formation of the clumps, eventually the self-gravity should become 
important, since the gravitational time scale is just a few times longer than the cooling 
time scale. As a result the initial clumps grow by both thermal and gravitational processes 
to become bound clouds. Hence, the thermal-gravitational instability should be considered 
(e.g., see Balbus 1986). 

In this paper we study the thermal-gravitational instability in a hot background whose 
physical parameters are relevant for protogalactic halo gas. Three-dimensional simulations 
were made. The simulations started with random Gaussian density perturbations of various 
power spectra. The evolution of gas under the influence of thermal-gravitational instability 
was followed up to the formation of gravitationally bound clouds. The physical properties 
of the clouds are examined in detail. Our goal is to study how self-gravity and gravitational 
interactions affect the physical and dynamical properties of clouds that condense initially via 
TI. Although our models would be yet too simple to be directly applied to the real situation, 
we try to extract the implications of our results on the formation of PGCCs in protogalactic 
halos. In §2 our models and numerical details are described. Simulation results are presented 
in §3, followed by summary and discussion in §4. 



2. Simulations 

2.1. Models for Protogalactic Halo 

A gas of Th = 1.7 x 10 6 K in a cubic, periodic simulation box was considered. This tem- 
perature corresponds to that of an isothermal sphere with circular velocity V c = 220 km s _1 , 
representing the halo of disk galaxies like the Milky Way. The fiducial value of the mean 
background density of hydrogen nuclei was chosen to be = 0.1 cm~ 3 . A case with higher 
density n h = 0.3 cur 3 was also considered to explore the effects of background density 
(Model D, see Table 1 below). For the primordial gas with an assumed ratio of He/H num- 
ber densities of 1/10, the gas mass density is given by ph = (2.34 x 10~ 24 g) n^. With 
Th = 1.7 x 10 6 K and — 0.1 cm -3 , the initial cooling time scale is t coo \ = 2 x 10 7 yrs. On 
the other hand, the free-fall time scale, or the gravitational time scale, is t grav = 1.4 x 10 s 
yrs, which is about seven times longer than the cooling time scale. Note that t coo i oc n^ 1 , 

— 1/2 

while tg rav oc n h . So cooling, compared to gravitational processes, becomes relatively more 
important at higher densities. The cooling length scale is given as / coo i = Ch ■ t coo \ = 4 kpc, 
where = 198 km s _1 is the sound speed. The simulation box was set to have the size 
L = 10 kpc = 2.5Z coo i- It was chosen to be large enough to accommodate a fair number 
of thermally unstable clouds of cooling length size, and so to get fair statistics of cloud 
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properties. 

To mimic density perturbations existed on a wide range of length scales inside the 
protogalactic halo, the initial density field was drawn from random Gaussian fluctuations 
with predefined density power spectrum. The density power spectrum was assumed to have 
the following power law 

P k = \5p k \ 2 4nk 2 oc k n , (1) 

where k is the three-dimensional wavenumber, k = [k 2 £ + k 2 + k 2 ) 1 / 2 . In order to explore how 
the initial density perturbations affect the formation and evolution of clouds, three types 
of density power spectrum were considered: white noise with n = as the representative 
case, as well as random fluctuation with n = 2 (Model R) and Kolmogorov spectrum with 
n = —5/3 (Model K). Only the powers with A;'s corresponding to wavelengths A < L/2 
were included. The amplitude of density power spectrum was fixed by the condition 5 rms = 
{Sp 2 ) 1 ^ 2 I (p) = 0.2. The initial temperature was set to be uniform, assuming isothermal 
density perturbations. The initial velocity was set to be zero everywhere in the simulation 
box. 



2.2. Numerical Details 

The gas-dynamical equations in the Cartesian coordinate system including self-gravity, 
radiative cooling and heating are written as 

^ + V-(pv) = 0, (2) 

d 

-(pv) + V • (pvv + pi) = -pV$, (3) 

dE 

— — h V • [(E + p)v] = — pv ■ V$ + T — A, (4) 
ot 

V 2 $ = 47rGp, (5) 

where E = (l/2)pv 2 +p/(7 — 1), A and T are the cooling and heating rates per unit volume, 
and the rest of the variables have their usual meanings. For the adiabatic index, 7 = 5/3 
was assumed. 

The hydrodynamic part was solved using an Eulerian hydrodynamics code based on the 
total variation diminishing (TVD) scheme (Ryu et al. 1993). The version parallelized with 
the Message-Passing Interface (MPI) library was used. 



The self-gravity, cooling and heating were treated after the hydrodynamic step. For the 
self-gravity, the gravitational potential was calculated by the usual Fast Fourier Transform 
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(FFT) technique. Then, the gravitational force was implemented in a way to ensure second- 
order accuracy as 

v n+1 = v n+1 '* - At n ■ V$ n+1/2 , (6) 

where v n+1 '* is the velocity updated in the hydrodynamic step and $ Tt+1 / 2 is the potential 
calculated with p n+1 / 2 = (p n + p n+l )/2. 

For the radiative cooling, the collisional ionization equilibrium (CIE) cooling rate for 
a zero-metalicity (primordial), optically thin gas in the temperature range of 10 4 K <T< 
10 7 - 5 K was adopted (Sutherland & Dopita 1993). The cooling function, L(T) = A/n 2 H , is 
plotted in Figure 1 as the solid line. With L(T) = for T < 10 4 K, it was assumed that the 
extra-galactic/stellar UV radiation photo-dissociates H 2 molecules in the protogalactic halo 
and so prohibits the gas from cooling below 10 4 K. Accordingly, the minimum temperature 
was set to be T min = 10 4 K. Note that the CIE cooling rate is higher than the cooling rate 
based on the non-equilibrium ionization (NEQ), which had been adopted in our previous 
two-dimensional simulations, especially near the H and He Lja line emission peaks (see 
Figure 1 of Baek, Kang & Ryu 2003). The higher cooling rate was chosen to accentuate the 
effects of cooling over those of gravitational processes. 

If H 2 molecules have formed efficiently via gas phase reactions enough to be self-shielded 
from the photo-dissociating UV radiation, or if the halo gas had been enriched by metals 
from first-generation supernovae, however, the gas would have cooled well below 10 4 K. In 
order to explore how the additional cooling below 10 4 K affects the formation and evolution 
of clouds, in a comparison model (Model C), the following mock cooling function in the range 
of 10 2 K < T < 10 4 K was adopted 

/ io 3 \ 

L H2 = 1.0 x 10~ 26 exp I -— J erg cm 3 s" 1 . (7) 

Although this mock cooling rate was designed to represent typical H 2 ro-vibrational line 
emissions for the gas with H 2 abundance of nn 2 /n H ~ 10~ 3 (Shapiro & Kang 1987), the 
exact amplitude and form of Lh 2 are not important in our discussion. Figure 1 shows Lh 2 
as the dashed line. The minimum temperature was set to be T min = 10 2 K in the C model. 

It was assumed that the background gas was initially under thermal balance and there 
existed a constant background heating equal to the cooling of the initial background gas, 
that is, 

r = L{T h )n\. (8) 

To prevent any spurious heating the highest temperature was set to be T max = T h . This ad- 
hoc heating was applied in order to maintain the temperature of the background gas with 
initial mean density at T ~ T h . This can be provided by several physical processes such as 
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turbulence, shocks, stellar winds and supernova explosions in protogalaxies. Without this 
heating the background gas would have cooled down in a few £ cool . 

It has been pointed out that the thermal conduction can affect profoundly TI in the 
ISM (e.g., McKee & Begelman 1990; Koyama & Inutsuka 2004). With the Spitzer value of 
thermal conductivity for ionized gas, k~5x 10~ 7 T 5 / 2 erg s -1 cm -1 K -1 (Spitzer 1979), the 
thermal conduction time scale is 

*- ~ 2 x 109 (it^m) (w^) (ii^) 2yrs - (9) 

It shows that t cond is much longer than t coo \ across the cooling length (Z cool = 4 kpc), so 
the thermal conduction is expected to be negligible. But the clouds finally emerged in our 
simulations have radius of R c ~ 150 — 200 pc (see §3.3). Even in such scale t con d is still a 
few times longer than, or at most comparable to, t coo i- Also those clouds are gravitationally 
bound, so the thermal evaporation should not be important. In addition, weak magnetic 
field, if it exists, may reduce significantly the value of thermal conductivity from the Spitzer 
value (e.g., Chandran & Cowley 1998), although we do not explicitly include magnetic field 
in this work. All together, it is expected that the thermal conduction does not play a major 
role in the regime we are interested in, and hence we ignored it in our simulations. 

2.3. Model Parameters 

Simulations were made with 1024 3 , 512 3 and 256 3 grid zones, allowing a uniform spatial 
resolution of Al = 9.8—39 pc. Simulations started at t — and lasted up to t cnd = 16— 20t coo \. 
This terminal time corresponds to ~ 2 — 3t grav , so gravitational bound clouds should have 
emerged by then. 

Total seven simulations are presented in this paper, which differ in numerical resolution, 
the power spectrum of initial density perturbations, cooling, and mean background density. 
Model parameters are summarized in Table 1. 

3. Results 

3.1. Evolution of Halo Gas 

We start to describe the results by looking at the global evolution of gas and its mor- 
phological distribution. Figure 2 shows the density power spectrum at different times in the 
S1024 model where oc constant initially. In the figure the dimensionless wavenumber is 
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given as k = L/A, which counts the number of waves with wavelength A inside the size L of 
simulation box. The power spectrum is presented in a way that 



With t coo i < tg rav , initially TI should work. The noticeable features in the early evolution 
of power spectrum are the followings: during t < 2t coo i the powers with k > 40 (or A < 
/ C ooi/16) are reduced significantly, and then the powers on those small scales grow back 
by ~ 4t coo i. The initial decrease of small scale powers is a consequence of initial isothermal 
density perturbations. The accompanying pressure fluctuations have generated sound waves, 
and those sound waves have ironed out the perturbations of small scales. The follow-up, 
fast growth of small scale powers is due to nonlinear behavior of TI. Although the linear 
growth rate is independent of scale (for A < / C ooi), the growth can be limited once the 
density increases and the cooling length becomes smaller than the perturbation scale. With 
the perturbation scale smaller than the cooling length, the further condensation progresses 
isochorically and the growth slows down. As a result, small scale clumps appear first. This 
point has been made previously, for example, by Burkert & Lin (2000). In addition, when 
perturbations get compressed by the background pressure, the density in the central region 
increases first and clouds form inside out. This contributes the fast growth of small scale 
powers too. By the end of this early TI stage, ~ 4t cool , the power spectrum peaks at 



After ~ 4t coo i, the self-gravity starts to play a role. Clouds grow through gravitational 
infall as well as compression by the background pressure. Eventually, massive clouds form 
through cloud-cloud collision, or merging among clouds (see below). During these stages, 
the powers grow over all scale. At the same time the peak shifts to smaller wavenumbers, 
reflecting the appearance of larger, massive clouds 

We note that with periodic boundary, once the power of the scale corresponding to the 
box size reaches nonlinear, the large scale clustering becomes saturated. The power spectrum 
in Figure 2 shows that in the S1024 model the scale which has gone nonlinear by the end is 
~ L/2 — L/3, indicating that the assumption of periodic boundary should not have affected 
our results significantly. In any case our major focus lies on the properties of individual 
clouds rather than their clustering. 

In order to show how clouds form and grow as well as how their distribution evolves, 
three-dimensional iso-density surfaces at four different times in the S1024 model are presented 
in Figure 3. As noted above, at the end of the early TI stage (~ 4t coo i) mostly small clouds 
appear. By the end of the follow-up stage of TI and gravitational growth, ~ 8t coo i, the clouds 
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A ~ ^cooi/50. 
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become larger. During late stage, the clouds become even larger but their number reduces, 
as a result of cloud-cloud mergers, which can be seen in the bottom two images. 

3.2. Identification of Clouds and Their Mass Function 

In order to study the properties of formed clouds they were identified using the algorithm 
CLUMPFIND described by Williams, De Geus & Blitz (1994). The algorithm basically tags 
cells around a density peak as the "cloud cells", if they satisfy the prescribed criteria of 
density and temperature. We chose the following criteria: p > 10po and T < 10 5 K. Here 
Po is the initial mean density. There is an arbitrariness in these threshold values. But the 
identification of clouds does not depend sensitively on the choices of threshold values, since 
clouds are well delineated by rather sharp jumps in density and temperature. In addition, 
we qualified only those with at least 3x3x3 cells or more as clouds. Once clouds were 
identified, their various quantities were calculated. 

The first row of Figure 4 shows the number of clouds, iV c , as a function of the cloud 
mass, M c , at the times same as those in Figure 3 in the S1024 model. By 4t coo i a significant 
number of clouds, N c ~ 3 x 10 3 , form through the growth of initial high density peaks by 
TI. The mass function at 4t cool is roughly Gaussian, since the initial density perturbations 
were drown from a random Gaussian distribution. During the follow-up stage of TI and 
gravitational growth, more peaks develop into clouds and at the same time they become 
massive. The mass function evolves roughly into the log-normal distribution, shown at 8t coo i- 
The log-normal distribution is a signature of nonlinear structure formation, as reported in 
various simulations (see, e.g., Vazquez-Semadeni 1994; Wada & Norman 2001). During late 
stage, more massive clouds develop through gravitational merging as pointed with Figure 3, 
and so the mass function extends to higher mass. The high mass tail of the mass function 
beyond the peak follows approximately a power-law distribution. When the mass function is 
fitted to dN/alM c oc M~ a for M c > 10 5 - 5 M , the value of a decreases from ~ 0.8 at 12t cool 
to ~ 0.6 at 16t coo i. 

3.3. Size, Density and Energetics of Clouds 

With the mass, M c , and volume, V c , of identified clouds, the effective radius was taken 
as R c = (3V C /An) x / 3 and the mean density was calculated as (p c ) = M c /V c . The second and 
third rows of Figure 4 show R c and (p c ) in the S1024 model. For energetics the kinetic energy 
relative to the cloud's center of mass, E K) the thermal energy, E T , and the gravitational 
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energy, E G) were calculated for identified clouds. The procedure by which the gravitational 
energy was calculated is presented in Appendix A along with a note of caution. The fourth 
row of Figure 4 shows the ratio of positive to negative energies, or the virial parameter 

— 2 • (Er + Et) 
\Eg\ 

The parameter f3 tells whether clouds are primarily pressure-bound (f3 > 2) or gravitationally 
bound (f3 < 2). Among the gravitationally bound clouds, the condition f3 ~ 1 indicates 
that they are approximately in virial equilibrium, although, strictly speaking, the condition 
applies only for a stable system where the external pressure is negligible and the moment of 
inertia does not change with time. 

During the early stage of TI (< 4t coo i), as small clumps grow isobarically, their density 
increases gradually, but (p c ) / p has not yet reached the isobaric factor of ~ (T h /T min ) ~ 100. 
The self-gravity is negligible with /?>2 during this stage. In the follow-up stage of TI and 
gravitational growth (< 8t coo i) , pressure-bound clouds develop fully. Those clouds have 
roughly (p c ) / p Q ~ 100, the maximum isobaric increase from the initial density, and the 
mean density is a bit higher in clouds with larger mass. Their radius scales roughly as 
R c ~ M°' 3 , which can be seen in the second row of Figure 4. The virial parameter has f3 > 2 
for all clouds, confirming that they are still gravitationally unbound. The pressure-bound 
clouds follow P oc M" 1 , which can be understood as follows. In those clouds, the thermal 
energy is dominant over the kinetic energy, i.e., E T ~ a few x E K . With T ~ T min in those 
clouds, the thermal energy scales as E T oc M c . On the other hand, the gravitational energy 
scales approximately as E G oc M^, since larger clouds have slightly more concentrated mass 
distribution. 

As clouds grow further primarily through merging in late stage, some of them become 
massive enough to be gravitationally bound. After that point clouds can be divided into 
two populations of distinct properties (see the last two times of Figure 4): pressure-bound 
clouds with smaller mass and gravitationally bound clouds with larger mass. The pressure- 
bound clouds have the properties similar to those found in the earlier stage. They have 
mean density (p c ) ~ 100 (pb g ), where (pb g ) is the mean background density. Note that the 
background density continues to decrease as more mass goes to clouds. By the end of the 
S1024 simulation, 16t coo i, only ~ 1/4 of gas mass remains in the background. So the mean 
density of the pressure-bound clouds decreases with time. Those pressure-bound clouds have 
(3 > 2 and follow (3 oc M c _1 . The gravitationally bound clouds appear first at ~ 9t coo i or 
~ 1.3t grav in the S1024 model. In those clouds the self-gravity enhances the density and the 
mean density of the clouds reaches up to 10 4 (pb g ) or even higher by 16t coo i . The fourth raw of 
Figure 4 shows that the gravitationally bound clouds are approximately in virial equilibrium 
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with f3 ~ 1. In addition, we found that in those clouds the kinetic energy is not smaller but 
sometimes larger than the thermal energy as expected in a virialized system, and both the 
gravitational energy, E G) and the positive energy, E K + E T , scale as M c 2 . 

Two points are noted on the gravitationally bound clouds. 1) Although these clouds 
have P ~ 1, they are not in steady-state. The clouds lose the positive energy rather quickly 
through cooling and contract further. But at the same time the cloud mass continues to 
grow through merging. 2) The gravitational energy of these gravitationally bound clouds 
scales as Eq oc M 2 , because their radius is in a relatively narrow range of 150 — 200 pc 
regardless of their mass, as shown in the second row of Figure 4. Since there is no physically 
obvious reason why the clouds with different mass should have similar radii, it should be 
understood as the result of dynamical evolution. However, their radius is somewhat larger 
at 16t coo i than at 12t cool . This is related to the increase of angular momentum in the clouds 
with time (see the next subsection). 

The distinction between the pressure-bound and gravitationally bound clouds in late 
stage can be understood with the critical mass 

M crit = l.l 8 f^) 2 G- 3 / 2 p- 1/2 , (12) 
\Hm H J 9 

which is the maximum stable mass for an isothermal sphere confined by the background 
pressure pb g (see e.g., McCrea 1957; Kang, Lake & Ryu 2000). For T c = 10 4 K and the 
background pressure four times smaller than the initial pressure (due to decrease in the 
background density), the critical mass is M crit ~6x 10 6 M , which coincides well with the 
transition mass scale in the last two times of Figure 4. Incidentally, this mass is similar to 
the Jeans mass of clouds, which is given as 

Mj = p\ 3 j = 6 x 10 6 ( ^-V /2 ( ^-^y V2 , (13) 

where Aj is the Jeans length (see e.g., Spitzer 1979). 

Although the simulated model would be too simple to represent a real protogalactic halo, 
it is tempting to regard the gravitationally bound clouds as possible candidates for PGCCs. 
That is, some of them may further cool down below 10 4 K either by UV self-shielding of 
H 2 molecules or by self-enrichment of metals due to first generation Type II supernovae, 
possibly leading to star formation. If about ten of them turn into globular clusters, their 
number density would be ~ 0.01 clusters/kpc 3 . However, the typical size of globular clusters 
~ 10 pc or so. So the clouds should contract further by more than a factor of 10. But the 
further collapse is controlled by the rotation of clouds (see §3.4 below). In addition, the 
typical mass of globular clusters is ~ 10 6 M Q or so. So if these gravitationally bound clouds 
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were to become globular clusters, the star formation efficiency should be ~ 10% or so with 
~ 90% of their mass dispersed back to protogalactic halo. 

It was pointed out by Truelove et al. (1997) that "artificial fragmentation" due to errors 
arising from discretization occurs in numerical simulations with self-gravity. They argued 
that the artificial fragmentation can be suppressed if resolution is maintained high enough 
that, for instance, the "Jeans number" Al/Xj < 0.25 or so for isothermal collapses. Here Xj 
is the local Jeans length. We found that although the constraint was not complied in a few 
high density zones, Al/Xj was kept to be always smaller than 0.4 in the S1024 simulation. In 
any case, we followed only up to the formation of bound clouds, not the subsequent evolution 
leading to fragmentation of those clouds. So no obvious fragmentation was observed. 



An important property of clouds that controls the dynamical state and affects the even- 
tual fate is their rotation. To quantify it, the angular momentum of clouds relative to their 
center of mass, J c , was calculated. The bottom row of Figure 4 shows the specific angular 
momentum of clouds, j c = J c /M c , in the S1024 model. 

In pressure-bound clouds, rotation plays a minor role in their dynamical evolution, 
since the rotational energy is smaller by an order of magnitude than the thermal energy. 
But we found that rotation can be dynamically important in gravitationally bound clouds. 
The specific angular momentum of those gravitationally bound clouds is larger than that 
of pressure-bound clouds, which can be seen clearly at 16t coo i . In fact the gravitationally 
bound clouds with a same M c have higher angular momentum at later time. For example, 
the clouds with M c ~ 10 7 - 5 M Q have a few times larger j c at 16t coo i than at 12£ coo i- This is 
because the clouds have grown mostly through merging in late stage and they have gained 
angular momentum through merging as well as tidal torque. Its direct consequence is that 
the clouds at later stage have larger radius, as noted in the previous subsection. We note 
that the clouds in our simulations are ever evolving. So it is not meaningful to define the 
canonical properties of clouds such as R c and function of M c . 

The rotation of gravitationally bound objects is often characterized by the dimensionless 
spin parameter (Peebles 1969) 



3.4. 



Rotation of Clouds 



X s 
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It measures the degree of rational support of the systems with 



centrifugal 



gravity 



Rl GM C 




(15) 



In the context of cosmological iV-body simulations, it has been shown that the typical value 
of \ s for gravitationally bound objects which acquired angular momentum through gravita- 
tional torque is ~ 0.05 (e.g., see Barnes & Efstathiou 1987). In Figure 5 red squares represent 
the spin parameter for the gravitationally bound clouds with M c > 10 7 M at 16t coo i in the 
S1024 model. It is clear that A s of our gravitationally bound clouds is significantly larger 
than 0.05. All except one have A s > 0.05 and the mean value is (A s ) ~ 0.3. It is because the 
clouds have gained angular momentum through merging as well as torque. 

With such large values of \ s , the degree of rotational support of the gravitationally 
bound clouds should be already substantial. Hence, rotation should be the key parameter 
that determines whether some of those clouds could become PGCCs and collapse further to 
globular clusters. For instance, in a dissipative collapse that conserves angular momentum, 
the clouds with A s ~ 0.3 can contract only by a factor of a few, before they become completely 
rotationally supported and disk-shaped. However, yet it is not clear whether we should 
conclude that those clouds with A s ~ 0.3 can not evolve into globular clusters. It is because 
we can not rule out the possibility that the clouds may be able to lose most, say ~ 99%, of 
their angular momentum, while they collapse and lose ~ 90% of their mass. For instance, 
in the so-called self- enrichment model of globular cluster formation, first generation Type 
II supernovae govern star formation and the removal of residual gas (e.g., see Parmentier 
et al. 1999; Shustov & Wiebe 2000). In this model the star formation may have occurred 
preferentially at the core, and most gas in the outskirt with large angular momentum may 
have been blown out. 

In addition, as noted above, the gravitationally bound clouds found at an earlier time 
have smaller angular momentum. So if the clouds were detached from background and 
started to collapse earlier, the angular momentum restriction would be somewhat less severe. 
Also the clouds emerged from different environments could have smaller angular momentum 
(for instance, see §3.8). In any case, subsequent evolution of the gravitationally bound 
clouds is beyond the scope of this paper, so the possible connection between these clouds 
and PGCCs should be left as a future study. 
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3.5. Shape of Clouds 

Shape is another property of clouds that reflects their dynamical state. In order to 
quantify it, we examined the shape parameters defined as 

b c 
q = - and s = -, (16) 
a a 

with each cloud fitted to triaxial ellipsoid with axes of a > b > c. The shape parameters 
have been commonly used to study clumps in numerical simulations (e.g., Curir, Diaferio 
& de Felice 1993; Gammie & Lin 2003). In order to find them, first the moment of inertia 
tensor, 

W^*' (17) 

was constructed for each cloud. Here x is the displacement relative to the cloud's center of 
mass and the integral is taken over the cloud volume. Then, from the three eigenvalues of 
the tensor, I n > 7 2 2 > ^33, the shape parameters were calculated as 

The clouds with s ~ q < 1 are of prolate shape and the clouds with s < 1 and q ~ 1 are of 
oblate shape, while triaxial clouds have s < q < 1. Spherical clouds have s ~ q ~ 1. 

In the previous study of isolated, thermally unstable clouds using two-dimensional sim- 
ulations in cylindrical geometry, we showed that the cloud shape changes in the course of 
the evolution (Baek, Kang & Ryu 2003). The degree of oblateness or prolateness is en- 
hanced during the initial cooling phase, as expected. But it can be reversed later due to the 
supersonic infall along the direction perpendicular to the initial flatness or elongation. 

In Figure 6 dots represent clouds in the q — s plane at four times in the S1024 model. 
Lines divide the domain into three regions of roughly prolate (left), triaxial (middle) and 
oblate (right) shapes. The panel for 4t coo i indicates that most clouds formed as a result of 
TI at early stage are preferentially of prolate shape. This is can be understood from the 
fact that filaments are the morphology dominant next to knots of clouds. By 8t coo i some of 
prolate clouds have been transformed to be oblate, as the result of the shape reversal, which 
was observed in the two-dimensional study. However, the figure shows that in late stage 
clouds tend to shift back to be prolate again. This is because gravitational merging results 
preferentially in clouds of elongated prolate shape. 

The gravitationally bound clouds with M c > 10 7 M , shown in Figure 5, are marked 
with red circles in the panel for 16t coo i- It is interesting to note that unlike most clouds, 
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those massive clouds are preferentially of oblate shape. It is because those clouds have large 
angular momentum and a substantial degree of rational support. 



Convergence is an important issue in any numerical simulations. We tested it by compar- 
ing the results of simulations with different resolution, i.e., S1024, S0512 and S0256 models 
(see Table 1). Figure 7 shows the density power spectrum of the S0512 and S0256 models, 
which can compared with that of the S1024 model in Figure 2. Note that in our simulations 
the amplitude of the initial power spectrum was set by the condition 



with 8 ims = 0.2 for all resolutions. Simulations with different resolution covers different 
range of wavenumbers. So the simulations of lower resolution started with larger amplitude, 
as shown in the figure. Otherwise, the evolution of the power spectrum looks similar in 
all three models. For instance, the two most noticeable features in Figure 2, i.e., the initial 
decrease and follow-up fast growth of small scale powers, are also present in Figure 7. But an 
interesting point is that the scale that suffered the initial decrease is insensitive to resolution, 
because it was caused by smoothing due traveling sound waves. On the other hand, the scale 
of the peak in the power spectrum after the follow-up fast growth does depend on resolution. 
As a matter of fact, the peak at 4t coo i in the S0256 model occurs at the scale almost four 
times larger than in the S1024 model. This indicates that the small scale growth in our 
simulations was limited by resolution, as expected. 

Figure 8 shows the mass function, virial parameter and specific angular momentum of 
clouds in the three models. The left panels show the results at early TI stage, while the 
right panels show the results at late merging stage. For early stage a different time was 
chosen for each model so that the density power spectrum has a similar amplitude. In late 
stage, however, the differences caused by the initial amplitude of power spectrum become 
insignificant, so 15 t coo i is chosen for all three models. We see that the mass function have 
been converged for massive clouds, i.e., those with M c > 10 6 M Q at the early stage and 
those with M c > 10 7 M Q at 15 t coo \. However, the number of smaller mass clouds depends 
on numerical resolution, as expected. With the minimum number of 3 3 zones for identified 
clouds, the minimum mass scales M min oc (A/) 3 . On the contrary, at the early stage the 
virial parameter is larger in lower resolutions for clouds of all mass. It is because the clouds 
are less compact in lower resolution, and so their gravitational energy is smaller. But as 
the clouds grow more massive and larger in late stage, the difference in the virial parameter 
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(19) 



- 15 - 



becomes smaller. The angular momentum of clouds is somewhat smaller in lower resolution 
at the early stage. It is partly because the time of the plot is different in different models. 
However, in late stage the angular momentum becomes comparable in all three models. 

Overall, the formation of small mass clouds, initially through TI and subsequently by 
compression and infall, was affected by resolution in our simulations. However, we found that 
the massive, gravitationally bound clouds, which have formed mostly through gravitational 
merging, have the properties which are almost converged. 

3.7. Effects of Initial Power Spectrum 

The effects of different initial perturbations on the formation and evolution of clouds 
were examined by comparing the results of simulations with different initial density power 
spectrum, i.e., the K0512 and R0512 models, to those of the S0512 model (see Table 1). 
The K0512 model started with more power on larger scales, while the R0512 model started 
with more power on smaller scales. Figure 9 shows the density power spectrum of the K0512 
and R0512 models. In both models, small scale powers suffered the initial decrease, as in 
the S0512 model. Especially most of the powers in k > 40 were erased substantially in the 
R0512 model. Hence, the overall growth was delayed in the R0512 model. On the other 
hand, the growth proceeded faster in the K0512 model, with more powers on large scales in 
the beginning. 

Figure 10 shows the total number of clouds as a function of time in the three models 
with different initial density power spectrum. The overall evolution is similar; the number 
of clouds increases during the TI and follow-up growth stages, but eventually decreases as 
gravitational merging progresses. But as noted above, the K0512 model evolves first and the 
S0512 and R0512 models follow. So clouds form from t ~ 2t coo i in K0512, from t ~ 5t coo i in 
S0512 and from t ~ 10t coo i in R0512. However, an interesting point is that the maximum 
number of clouds is about the same with iV max ~ 3000 in all three models. 

In Figure 11 three-dimensional iso-density surfaces are plotted at two sets of times 
in the three models. Different times were chosen in different models, since the formation 
and evolution of clouds proceeds differently. The upper panels show the surfaces when the 
number of clouds is highest, i.e., at 3t coo i for K0512, 6t coo i for S0512 and llt coo i f° r R0512. 
The lower panels show the surfaces after the number of clouds have decreased a little bit due 
to merging, i.e., at 4t coo i for K0512, at 10t coo i for S0512 and at 19t coo i for R0512. The most 
noticeable feature is that the distribution is more "filamentary" in the K0512 model, but less 
"filamentary" in the R0512 model, than in the S0512 model. It is because that the initial 
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large scale powers, which were largest in the K0512 model but almost absent in the R0512 
model, have been developed into significant coherent structures in the cloud distribution. 

Figure 12 shows the mass function, virial parameter and specific angular momentum of 
clouds in the three models. The same two sets of times as those in Figure 11 were chosen. 
Clouds in the K0512 (R0512) model are slightly more (less) massive in the left panels and 
slightly less (more) massive in the right panels. However, considering the difference in 
the plotted time, the cloud mass function should be regarded as reasonably similar in the 
three models. On the other hand, the virial parameter of pressure-bound clouds follows the 
same diagonal strip in all three models. There is a spread in the distribution of angular 
momentum, again partly because the plotted time is different in different models. But the 
angular momentum of the clouds in high mass tail is similar. So we conclude that the 
properties of individual clouds, especially for massive clouds, are not sensitive to the initial 
perturbations, while the spatial distribution of clouds reflects the slope of initial density 
power spectrum. 



3.8. Effects of Different Density and Cooling 

The effects of gas density on the formation and evolution of clouds were examined with 
the D0512 model, which has the background density 3 times larger than that of the S0512 
mode. Otherwise the two models are identical (see Table 1). With t coo \/t gra , v oc n^ 2 , 
the gravity is relatively less important in the D0512 model. Figure 13 compares the mass 
function, virial parameter and specific angular momentum of clouds in the D0512 model 
(red lines and dotes) with those in the S0512 model at an early stage of TI (5t coo i) an d at 
a late stage of merging (15t coo i), respectively. In the D0512 model there are more clouds 
with larger mass, reflecting the higher background density. Yet, the virial parameter is 
almost identical in the two models, indicating that the background density is not important 
in determining the energetics of individual clouds. As a matter of fact, we found that the 
properties of clouds, except the angular momentum, are not sensitive to the background 
density. The angular momentum of the clouds identified in the early stage of TI is similar 
in the two models. But in late stage the angular momentum is noticeably smaller in the 
higher background density model. The same trend is also obvious in the spin parameter 
of gravitationally bound clouds, which is shown for M c > 10 7 M in Figure 5. While the 
spin parameter for the S0512 model (blue triangles) does not differ much from that for the 
higher resolution S1024 model (red squares), the spin parameter for the D0512 model (green 
circles) is substantially smaller with the median value of A me d ~ 0.12. This is because the 
clouds formed in higher background density have experienced relatively less tidal torque and 
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gravitational merging, through which they have acquired angular momentum. 

The effects of possible cooling further below 10 4 K were examined with the C0512 model, 
which includes a mock H 2 cooling given in equation (7). Otherwise it is same as the S0512 
model (see Table 1). Figure 13 compares the mass function, virial parameter and specific 
angular momentum of clouds in the C0512 model (black lines and dotes) with those in the 
S0512 model. The cloud properties are similar overall, except the smaller thermal energy 
in the C0512 model, which is expected from the additional cooling. Because the thermal 
energy counts for most of the positive energy, especially in pressure-bound clouds, the virial 
parameter is smaller in the C0512 model. But in massive, gravitational bound clouds, the 
kinetic energy is comparable to or sometimes larger than the thermal energy, as noted in 
§3.3. So the effects of the additional cooling is less important in those clouds. 

To quantify how much the additional cooling changes the thermal state of gas, we 
compare the mass distribution, f(T), for the S0512 and C0512 models in Figure 14. The gas 
with T ~ 10 4 K in S0512 spreads below 10 4 K in C0512, as expected. In addition some of 
the gas with 10 4 <T<2x 10 4 K has cooled below 10 4 K. But the additional cooling does 
not affect much the gas of higher temperature. Even with this additional cooling the mass 
fraction peaks still at T ~ 1.5 x 10 4 K and most of the "cloud gas" has 10 3 < T < 10 5 K, 
not only because the assumed mock H2 cooling is inefficient, but also because some of the 
gas has been reheated by shocks and compression. 

Here we should note that with lower temperature, the clouds in the C0512 model have 
smaller Jeans mass and could go through further fragmentation. But with a fixed grid 
resolution, the simulation could not follow it. In fact, in the C0512 model, the Jeans number, 
Al/Xj, reached up to > 1 in the center of some gravitationally bound clouds. However, as 
noted above, we did not intend to follow such fragmentation in this study. 

4. Summary and Discussion 

We study the role of self-gravity and gravitational interactions in the formation of clouds 
via thermal instability (TI) through three-dimensional hydrodynamic simulations. We con- 
sidered the gas in protogalactic halo environment with T ~ 1.7 x 10 6 K and n ~ 0.1 cm~ 3 in a 
periodic box of 10 kpc and followed its evolution for up to 20 cooling time. We adopted ide- 
alized models in which a static, non- magnetized gas cools radiatively with initial isothermal 
density perturbations. A radiative cooling rate in ionization equilibrium for an optically thin 
gas with the primordial composition was used. In addition, an ad hot heating was included, 
which emulates feedbacks from stellar winds, supernovae, turbulence and shocks in order to 
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maintain the thermal balance of background gas. 

The main results can be summarized as follows: 1) Clouds form first on scales much 
smaller than the cooling length as a result of non-linear behavior of TI. 2) Those small 
clouds grow through compression by background pressure as well as gravitational infall, 
but eventually they merge by gravity to become gravitationally bound objects. 3) The 
gravitationally bound clouds have acquired angular momentum through merging as well as 
tidal torque. So they have high angular momentum with the spin parameter of (A s ) ~ 0.3 
or so. 4) The spatial distribution of clouds depends on initial perturbations, for instance, 
the slope of initial density power spectrum, but the properties of individual clouds are not 
sensitive to that. 

We note that the realistic picture of thermal-gravitational instability that in protogalac- 
tic halos should be more complex than in the numerical models considered here. Some of 
key aspects include: 1) The gas in real protogalactic halos is likely in a chaotic state induced 
during the formation of halos themselves. The chaotic flow motions would have suppressed 
the early formation of clouds via TI, but increased collisions of clouds once formed. 2) Pro- 
togalactic halos have their own structures, but the effects of those structures were ignored. 
For instance, the tidal torque exerted by the halo and/or the rotation in disk would have 
suppressed the formation of clouds. 3) There are emerging evidences that magnetic field ex- 
isted even in the early galaxies where the oldest stars formed (see, e.g., Zweibel 2003). Then, 
undoubtedly the magnetic field should have affected the formation and properties of clouds 
profoundly. 4) It is well known that when a hot gas cools from T > 10 6 K, it recombines out 
of ionization equilibrium because the cooling time scale is shorter than the recombination 
time scale (Shapiro & Kang 1987). However, details of the cooling such as non-equilibrium 
ionization and H 2 and metal cooling below 10 4 K could have only minor effects on our main 
conclusions. 

Although the numerical models are rather idealized to facilitate simulations, our results 
should provide crude insights on the formation of PGCCs in protogalactic halos. The gravita- 
tionally bound clouds in our simulations have mass M c > 10 7 M Q and radius R c ~ 150 — 200 
pc. If some of them evolved into PGCCs and became globular clusters, they should have 
lost ~ 90% of their mass with ~ 10% of star formation efficiency, and at the same time 
they should have collapsed by a factor 10 or so. But the further collapse would not have 
been straightforward because of large angular momentum, unless their angular momentum 
was removed very efficiently along with mass during the star formation phase. Such removal 
of angular momentum may not be impossible, but following it is beyond the scope of this 
numerical study. It should be studied with simulations that have resolution high enough 
to follow the fragmentation of clouds and the subsequent formation of stars inside PGCCs. 
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However, we make the following note of caution. The clouds in our simulations continue 
to grow mostly through merging in late stage, and thus their mass and angular momentum 
increase in time. So any simulations of isolated clouds to study the ensuing evolution of 
PGCCs could be misleading. 

We should point that there is a caveat in our argument for the large angular momentum 
of gravitationally bound clouds. The angular momentum was acquired mostly through grav- 
itational processes, i.e., tidal torque and merging. So if the clouds formed in an environment 
where the gravitational processes are less important, they would have acquired less angular 
momentum. For instance, we showed that the clouds formed in a denser background have 
smaller angular momentum. Hence, if protogalactic halos consisted of smaller halo-lets and 
clouds formed in shocked regions after collisions of halo-lets, as suggested by e.g., Gunn 
(1980), they could have smaller angular momentum. 

We thank the anonymous referee for constructive comments. This work by CHB, HK 
and JK was supported by KOSEF through Astrophysical Research Center for the Structure 
and Evolution of Cosmos (ARCSEC). This work by DR was supported by Korea Research 
Foundation Grant (KRF-2004-015-C00213). Numerical simulations were performed using 
"Linux Cluster for Astronomical Calculation" of KASI- ARCSEC. 



A. Corrected Potential 

In our simulations, the gravitational potential for the gas-dynamical equation in §2.2 
was calculated by using the FFT method, <3> F ft(j)- Then the gravitational force could be 
correctly calculated by differentiating this potential on the grid. However the use of $fft(j) 
in calculating the gravitational energy of clouds, Eg, in §3.3 ends up a large error, because 
$fft(^) 1S undetermined by an integral constant. In principle, the gravitational potential of 
isolated clouds can be precisely calculated by the direct double integration over cloud volume. 
However, the computational cost of this method is prohibitively expensive, especially for 
gravitationally bound clouds in the 1024 3 simulation, since they occupy typically ~ 10 4 or so 
grid zones. One the other hand, the direct integration does not take account of contributions 
from the rest of mass in the simulation box as well as the periodic mass distribution. But 
we found that those contributions are small, especially for massive, gravitationally bound 
clouds. 

As an effort to estimate the gravitational energy of clouds more accurately, we devised 
a method which calculates and uses the corrected potential as follows: 

1) The position, r min , where $fft(^*) has the minimum value, is found for each cloud, and 
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then the potential at the position is calculated by the direct integration, 

$direct(r min ) = —G I -J^-L-dV'. (Al) 

J y 1 mm | 

2) The difference between ^direct and $fft at r min is calculated for each cloud, 

C = $direct(r mi n) - $FFT(r min ). (A2) 

3) Then the corrected potential for each cloud is calculated by 

<Wr) = <Wf) + a (A3) 

4) Finally the gravitation energy of each cloud is calculated as 

E G = J l -p{r)<S> colI {r)dV. (A4) 

Figure 15 demonstrates the motivation of our effort. Here is the gravitational energy 
of clouds calculated by the direct double integration, while ^fft and -E C orr were calculated 
using $ FFT and $ CO rr, respectively. The errors for the S0256 model are shown, where the 
double integration could be done with a reasonable computation time. E corr tends to agree 
with E N 2 better than _E F ft- With _E corr the error is within ~ 20% or so for a substantial 
fraction of clouds. However, the estimation of gravitational energy could be easily off by a 
factor two or even larger with E FFT . We used -E corr for the gravitational energy in §3.3. 
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Table 1. Model Parameters for Simulations 1 



Model 


No. of grid zones t cnd (t coo i) b W (£ g mv) c 


T 1 (~W\ 


Pk n h ( 


[cm 


•3\d 
) 


S1024 


1024 3 


16 2.29 


10 4 


const 


0.1 




S0512 


512 3 


20 2.86 


10 4 


const 


0.1 




S0256 


256 3 


20 2.86 


10 4 


const 


0.1 




K0512 


512 3 


20 2.86 


10 4 


oc k~'i 


0.1 




R0512 


512 3 


20 2.86 


10 4 


oc k 2 


0.1 




C0512 


512 3 


20 2.86 


10 2 


const 


0.1 




D0512 


512 3 


20 1.65 


10 4 


const 


0.3 




Simulation box size L 


= 10 kpc in all models. 










W 

''cool 


= 2 x 10 7 yrs in 


the models with rih = 0.1 cm" 


" 3 , and t cool 


= 6.7 x 10 6 


yrs 


in the 


model with n h = 0.3 cm 


3 










c t 

''grav 


= 1.4 x 10 8 yrs in the models with nh — 0.1 cm 




= 8.1 x 10 7 


yrs 


in the 


model with n h = 0.3 cm 


3 










Ph - 


(2.34 x 10- 24 g)n h with n(He)/ra(H) = 0.1. 











-24- 




Fig. I. — Cooling function, L(T) = A/n 2 H , of the collisional ionization equilibrium (CIE) 
cooling for a gas with zero-metalicity (solid line). The dotted line is for the mock H 2 cooling 
adopted in Model C. 
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Fig. 2. — Evolution of the density power spectrum in the S1024 model. Dots represent the 
initial power spectrum at t — and lines are the power spectrum at t coo \, 2 t coo \, 3 t coo i, • • •, 

16 tcool- 
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Fig. 3. — Iso-density surfaces inside the full simulation box of size 10 kpc at 4 £ coo i (top- 
left), 8 t coo i (top-right) , 12 £ cool (bottom-left) , and 16 £ coo i (bottom-right) in the S1024 model. 
Green surfaces corresponds to 10p , yellow surfaces corresponds to 10 2 p and red surfaces 
corresponds to 10 3 po- Here po is the mean initial density. 
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Fig. 4. — Differential number of clouds, dN c /d(\ogM c ), effective radius, R c , mean density, 
(p) c , energy ratio, 2(E T + E K )/\E G \, and specific angular momentum, j c , as a function of 
cloud mass, M c , at four different times in the S1024 model. Here, j c is in the cgs units 
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Fig. 5.— Spin parameter, A s = J\Eq /2 \/GM 5 / 2 , as a function of cloud mass, M c , for clouds 
with M c > 10 7 M . Red squares are for the S1024 model, blue triangles are for the S0512 
model, and green circles are for the D0512 model at 16 t C ooi- 
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Fig. 6. — Shape parameters, q = b/a and s = c/a, at four different times in the S1024 model. 
The principal axes of clouds are defined such that a > b > c. The domain is divided into 
three regions containing the clouds of prolate shape (s ~ q < 1, left), and the clouds of 
oblate shape (s < 1 and q ~ 1, right) and the clouds of triaxial shape (s < q < 1, middle). 
Red circles represent the clouds with M c > 10 7 M Q . 
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Fig. 7. — Evolution of the density power spectrum in the S0256 and S0512 models. Dots 
represent the initial power spectrum at t — and lines are the power spectrum at t coo \, 2 t coo \, 
3 ^cooh • • •) 20 t coo \. 
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Fig. 8. — Differential number of clouds, dN c /d(\ogM c ), energy ratio, 2(Er + Ek)/\Eq\, and 
specific angular momentum, j c , as a function of cloud mass, M c , in the three models with 
different resolution. Left panels show the quantities at 5 t coo i, 6 t coo i; an d 7 t coo i for the S0256, 
S0512, and S1024 models, respectively, and right panels show the quantities at 15 t coo i for 
all three models. Here, j c is in the cgs units. 
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Fig. 9. — Evolution of the density power spectrum in the K0512 and R0512 models. Dots 
represent the initial power spectrum at t — and lines are the power spectrum at t coo \, 2 t coo \, 
3 ^cooh • • •) 20 t coo \. 
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Fig. 10. — Time evolution of the number of clouds, N c , in the three models with different 
initial density power spectrum. 



Fig. 11. — Iso-density surfaces inside the full simulation box of size 10 kpc in the K0512 
model (left images), in the S0512 model (middle images), and in the R0512 model (right 
images). Top images are at 3 t coo i in the K0512 model, at 6 t coo i in the S0512 model, and at 
11 t CO oi in R0512 model. Lower images are at 4 t coo i in the K0512 model, at 10 t coo i in the 
S0512 model, and at 19 t coo i in the R0512 model. Green surfaces corresponds to 10po, yellow 
surfaces corresponds to 10 2 p and red surfaces corresponds to 10 3 p , same as in Figure 3. 
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Fig. 12. — Differential number of clouds, dN c / d(\og M c ) , energy ratio, 2(Et + Ek)/\Eg\, and 
specific angular momentum, j c , as a function of cloud mass, M c , in the three models with 
different initial density power spectrum. Left panels are at 3 t coo i, 6 t coo i, and 11 t coo i for the 
K0512, S0512, and R0512 models, respectively, and right panels are at 4 t cool , 10 t coo i, and 
19 t coo i for the K0512, S0512, and R0512 models, respectively. Here, j c is in the cgs units. 



-36- 



O 

o 

o 



o 



H 
X 

cv 

QjO 
C 



C 




6 



8 



'M 



8 

log(M c 

Fig. 13. — Differential number of clouds, dN c /d(logM c ), energy ratio, 2(Et + Ek)/\Eq\, and 
specific angular momentum, j c , as a function of cloud mass, M c , in the C0512 and D0512 
models as well as in the S0512 model for comparison. Left panels are at 5 t C oob and right 
panels are at 15 t coo i, respectively, for all three models. Here, j c is in the cgs units. 
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Fig. 14. — Mass fraction, f(T), of gas as a function of temperature at three different times 
in the S0512 and C0512 models. 
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Fig. 15. — Error in the estimation of gravitational energy by the FFT potential (bottom 
panels) and the corrected potential (top panels) at 10 t coo i (left panels) and 15 t coo i (right 
panels) in the S0256 model. See the Appendix A for details. 



